Casimir-Polder energy level shifts of an out-of-equilibrium particle near a microsphere 
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Rydberg atoms and beams of ultracold polar molecules have become highly useful experimental 
tools in recent years. There is therefore a need for accessible calculations of interaction potentials 
between such particles and nearby surfaces and structures, bearing in mind that the particles are 
far out of thermal equilibrium with their environment and that their interaction is predominantly 
non-retarded. Based on a new perturbative expansion with respect to the inverse speed of light and 
the inverse conductivity, we derive a simple, closed- form expression for the interaction potential (i.e., 
the particle energy level shifts) of a particle and a metallic sphere that is is accurate at better than 
1% level for typical experimental set-ups at room temperature and above, and off by no more than 
a few percent at any temperature including zero. Our result illuminates the influence of retardation 
and imperfect conductivity and the interplay of these effects with geometry. The method developed 
for the present study may be applied to other, more complex geometries. 

PACS numbers: 31.30.jh, 12.20.-m, 34.35.-|-a, 42.50.Nn 



I. INTRODUCTION 

Recent times have witnessed a blossoming of experi- 
mental set-ups in which the detailed interaction of par- 
ticles with nearby surfaces, in particular the Casimir- 
Polder (retarded van der Waals) [l| interaction, is im- 
portant. Some such systems are typically far out of ther- 
mal equilibrium, such as Bose-Einstein condensates in 
magnetic traps close to surfaces, beams of cold polar 
molecules and Rydberg atoms. We shall focus on the 
two latter categories herein. For example, the interaction 
between Rydberg atoms '^ and surfaces are essential to 
the understanding of the behavior of Rydberg atoms in 
vapour cells [3| and near atom chips [4], systems which 
have already been investigated in several experiments. 
Various suggested mechanisms for quantum information 
processing also involve Rydberg atoms close to metallic 
surfaces 0-0|- Moreover, beams of cold polar molecules 
have already been put to use in a range of experimen- 
tal applications as reviewed in Refs |8l4lOl|. For instance, 
trapping of cold CO molecules near atom chips using elec- 
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trie traps has recently been realized 

We derive herein a simple closed form expression for 
the CP interaction between a particle and a sphere valid 
both for Rydberg atoms and cold molecules. Second only 
to the plane surface, the spherical geometry is arguably 
the most generically useful to consider in all microscopic 
applications. The microsphere is the standard vehicle in 
the rapidly progressing field of micromanipulation and 
photonics, a field closely bordering on atomic physics 
where CP forces are of importance. Using laser beams, 
microspheres can be trapped and pushed [12] and per- 
haps even pulled [13| for detailed manipulation, and are 
readily transported along optical fibres via the evanescent 
field \1^ . Microsphere optical resonators with extremely 



high Q-factors have been built, and are useful e.g. for 
low threshold lasing [l^. Our closed form expression, 
not involving the typical lengthy sums of Mie scattering 
coefficients, is immediately useful for direct insertion into 
numerical simulations of microsystems, as well as analy- 
sis of experimental data. 

Rydberg atoms and cold (ground state) polar 
molecules share two traits that set them sharply apart 
from the ground state or thermalized atoms which have 
typically been considered in the van der Waals and 
Casimir-Polder (CP) literature. Firstly, they are both 
far out of thermal equilibrium with their thermal environ- 
ment. Rydberg atoms have been excited to a high prin- 
cipal quantum number, far from an atom's thermalized 
state, which is almost identical to its ground state since 
excitation energies are large compared to k^T (room tem- 
perature assumed). The excitation energies of rovibria- 
tional states of molecules, in contrast, are small compared 
to /cbT, so a thermalized molecule significantly occupies 
a number of its energy eigenstates. Thus also ground 
state molecules at room temperature are far from thermal 
equilibrium. Systems with magnetic transitions exhibit 
similar properties [16[. 

Secondly, and for the same reason, the Casimir-Polder 
interaction between these particles and nearby surfaces is 
predominantly non-retarded, and retardation corrections 
due to the finite speed of light enter only as a correction. 
Typically, the retarded interaction regime stretches for 
tens and hundreds of micrometers for polar molecules 
and Rydberg atoms, respectively [l3,|la|: thus including 
the separations normally encountered in experiments and 
applications. 

Whereas thermal non-equilibrium initially complicates 
theoretical treatment, the non-retardedness of the inter- 
action introduces a significant simplification, allowing al- 



most surprisingly simple results to be achieved. The gen- 
eral theory for Casimir-Polder interactions of a particle 
in an arbitrary superposition of eigenstates was recently 
derived by some of us JlOf], and has since been applied 
to p lanar geometries for Rydberg atoms and molecules 
[ISL 120] • Strikingly, it was found that for non- retarded 
interaction with a flat metallic surface the interaction 
potential is virtually independent of temperature [2l|, 
a result that can be extended to arbitrary geometries 
|23 |. For a full theoretical background of different non- 
equilibrium CP scenarios the reader may additionally 
refer to Refs. |23l - l29| . The theories for different non- 
equilibrium situations, albeit apparently disparate, may 
be shown to concord as they should [301 ■ 

We consider the general situation of a par ticle whose 
eigenstates are |n). It was shown in Ref. [l^ that for a 
particle prepared in an arbitrary superposition of eigen- 
states \4>) = ^^Pn\n) with occupation probabilities p„, 
the Casimir-Polder potential may be written as a sum 
over transitions between pairs of eigenstates according 
to 



with 



U^ = y^^PnUn 



Un — 2_^ ^nk 



(1.1) 



(1.2) 



where the sum runs over all other eigenstates \k) to which 
there is an allowed dipole transition. 

Crucial to the understanding of the CP interactions of 
both cold molecules [2^ and Rydberg atoms 18] is the 
realization that only a few transitions turn out to give sig- 
nificant contributions. To wit, the important transitions 
were found to be those corresponding to the smallest dif- 
ference in eigenenergy AEf^n — E^ — E^ , i.e. the smallest 
transition frequency u^n = A£'fc„/?i or correspondingly 
the longest transition wavelength Afe„ = 27rc/a;fc„. For 
example, a Rydberg atom near a half-space prepared in 
an s-state of principal quantum number n, obtains sig- 
nificant contributions from transitions to the few differ- 
ent p-states of principal quantum numbers n and n ~ 1 
|18| : for ground state LiH molecules the only significant 
transition was to the lowest rotational state, whereas for 
YbF also the first vibrational state was required [20| . 
Transitions with larger AEkn could be ignored to a good 
approximation. 



Because of this fact the typical values of \kn for cold 
polar molecules and Rydberg atoms alike are usually 
much larger than the typical distance z from the par- 
ticle to a nearby body in experiments involving surfaces. 
In other words, for the dominating transitions \n) — >■ |fc). 
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hence the interaction is essentially non-retarded. 

Recently we found that in the non-retarded regime 
the CP interaction near a metallic half-space is virtu- 
ally temperature independent [21j. The thermal CP po- 
tential then agrees with its zero-temperature counter- 
part for all temperatures. This was later shown to be a 
reasonable approximation for bodies of arbitrary shape 
[22| . Temperature-dependent corrections to the zero- 
temperature potential were identified to stem from re- 
tardation and imperfect conductivity. The magnitude of 
the latter corrections were found to strongly depend on 
the body shape and curvature, demonstrating that ge- 
ometry and temperature are closely intertwined [31|. It 
is therefore necessary in practice to study different ge- 
ometries individually. 

The case of an atom interacting with a metal sphere 
to be studied is a prototype of a body with a curved 
surface. Various embodiments of the particle-sphere in- 
teractions at zero tem per ature have been treated by a 
number of authors |32h39| . In most of these works, the 
CP potential is obtained from a numerical computation 
which suffers from poor convergence at small curvatures. 
In contrast, we will derive an approximate analytical re- 
sult based on a perturbative expansion that is readily 
accessible while illuminating the impact of retardation 
and imperfect refiection as well as the interplay of these 
factors with geometry: We calculate the CP potential 
contribution from transition \n) -^ \k) under the rele- 
vant assumption of non-retarded interaction, Eq. (11.31) . 
which we quantify by a retardation parameter 
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<1 



(1.4) 



where r is distance from the particle to the sphere's cen- 
ter. We quote here the final result, to be derived below, 
for the thermal Casimir-Polder potential of an atom at 
distance r from the center of a metal sphere of radius R: 
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where we have introduced the dimensionless geometry 
parameter 



(j) = R/r. 



(1.6) 



The expression (|1.5|) is remarkably simple compared to a 
numerical evaluation of the starting equations, involv- 
ing infinite sums over Mie scattering coefficients. It 
holds when a; ^ 1, but Im{-y/e-}x0 — hn{y/e}ujknR/c 
still significantly exceeds unity, which is the case for 
good conductors in combination with typical values of 
X for the systems under consideration. The dots indi- 
cate higher-order contributions in the small parameters 
X and l/{xy/e). For definitions of the various quantities 
in Eq. (jl.Sp . see Sec. HIl below. 

In the following we derive the CP potential for a parti- 
cle near a metallic sphere including the leading correction 
for small x (retardation correction) and l/{xy^) (imper- 
fect reflectivity correction) , starting from the general CP 
theory for particles out of thermal equilibrium, which is 
outlined in Section |IT1 Our method may in principle be 
employed for any geometry to derive perturbative tem- 
perature corrections such as that presented herein. Ex- 
plicit corrections for the particle-sphere configuration are 
derived in Section Hill and analysed in Section PlI CI 



II. GENERAL FORMALISM 

The general formalism for the temperature-dependent 
CP force on a particle in an energy eigenstate is found 
in Ref. [131 . Here we shall restrict our attention to the 
special case of an isotropic particle. As explained above, 
the potential on a Rydberg atom or cold molecule alike 
takes form of a sum over just a few contributions from 
pairs of eigenstates which all have transition wavelengths 
in the same order of magnitude. It is sufficient therefore 
to consider a single transition \n) — >■ \k). Our consider- 
ations for one such transition will therefore hold for all 
relevant transitions, and all that is required in order to 
return to the full description of these particles is to sum 
the final result over the relevant transitions according to 
Eq. (|1.2I) . The generalisation to anisotropic particles is 
straightforward (cf. Refs [Ulil). 

This section reviews the general framework for treat- 
ing the small temperature correction in the non-resonant 
regime for arbitrary geometries, which w e ap ply below to 
a metal sphere. In the notation of Ref. [22], the CP po- 
tential of an isotropic non- magnetic particle in state \n) 
due to a possible transition to state |fc) splits naturally 
into a non-resonant (nr) and a resonant (r) part (l9| 



where the two parts are given by 






2fcBT|d„fc|2u;fc„ ^, riC,(r) 
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n(wA;n)R.er„^„(r). 



(2.1) 

(2.2a) 
(2.2b) 



Here, dfc„ — (fc|d|n) is the transition dipole matrix ele- 
ment and 
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(2.3) 



are the Matsubara frequencies. The photon number at 
frequency lo and temperature T is given by the Bose- 
Einstein distribution, 



n{(jj) 



1 



eyi^{fiw / k^T) 



-[n(-c^) + l]. (2.4) 



The function 



r.(r) 



lim trG(^)(r,r', 



(2.5) 



is given in terms of the scattering part G'-^' of the total 
dyadic Green's function satisfying 



V X V X 



— £(r,a; 



G(r,r',a;) = (5(r-r')l. (2.6) 



(1: unit tensor). The relative permittivity e(r,w) of the 
present bodies is isotropic and we have assumed the bod- 
ies to be non-magnetic. Due to causality, Fi^^ is real, 
being a generalized susceptibility evaluated at imaginary 
frequency, so in particular Fjo = Fp is real. 

In the following, we consider a single transition \n) — > 
|fc) and simplify our notation according to 

|dfc„P -^ |dp; uJkr, -> w; Unk{r) ^ U{r). 

Note that ui can be either positive or negative depending 
on whether the transition is upwards or downwards. 

Knowing that the potential is largely temperature 
independent through the non-retarded region, con- 
sider for now the regime in which the linear T- 
corrections becomes important, i.e., the spectroscopic 
high-temperature regime, 

/cbT ^ hu!. 

Here, the contribution of the lowest Matsubara frequency 
{j — 0) dominates in Eq. (I2.2a[) . so 
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The photon number in this regime is 

so the resonant potential (|2.2bp reads 
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Combining these results, we find for the full potential in 
the spectroscopic high-temperature regime that 
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(2.10) 



where AF^ = r^ - Tq. 

For comparison, in the zero-temperature hmit, in 
which the Matsubara sum becomes an integral accord- 
ing to standard procedures (e.g. the Euler-Maclaurin for- 
mula), one obtains 



C/(r) 



T=0 
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trrW(r) 
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(2.11) 



with Q{x) denoting the unit step function. 

In the nonretarded and perfect-conductor limits, we 
have ri^(r) ~ Rertj(r) ~ ro(r) [21] which implies 
ReAr^(r) ~ 0. In this case, both Eqs. ((2lB and (dHj) 
reduce to 
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(2.12) 



and the CP potential is independent of temperature 
throughout. 

In this article, we are interested in the corrections to 
the temperature-independent result (|2.12p that arise due 
to small violations of the non-retarded limit and perfect 
reflectivity. As seen from Eq. (|2.10p . they are governed 
by ReArt^(r). When all present macroscopic bodies are 
perfectly conducting (PC), then r,^ satisfies [24I 



^PC 



so that 
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0, (2.13) 



(2.14) 



is quadratic in w. This correction accounts for the fact 
that electromagnetic interactions are transmitted at the 
finite speed of light; we will refer to it as the retardation 
correction in the following. 

For an imperfect conductor, the corrections to the 
Green's function for small frequencies includes a second 
correction due to the frequency-dependence of the reflec- 
tivity of the bodies. We write 



AF„ 



AP"; 



ap: 



rofl. 



(2.15) 



due to retardation and reflectivity, respectively. When 
treating uj and e(w) as independent variables, the re- 
tardation correction AP™*- is the leading-order term in 
l/^e{uj) and next-to-leading in w; whereas the reflec- 
tivity correction APJ^°**' is the contribution sub-leading 
in l/y^e{uj) and leading in w. Note that the perfect- 
conductor limit |e(w)| — >■ c» does not commute with the 
nonretarded limit w ^> in this case. The incompati- 
bility of the two limits was first pointed out in Ref. [40| 
and is at the heart of the debate over the temperature 
correction to the Casimir effect [41| . For a metal body at 
typical frequencies and distances, the perfect-conductor 




FIG. 1: The geometry considered: a quantum particle pre- 
pared in eigenstate \n) outside a gold sphere. 



limit has to be performed before the nonretarded limit, 
see Sect, mil below. 

With the leading corrections to the Green's func- 
tion being given by Eq. (|2.15p . the thermal CP poten- 
tial (|2.12l) can be given as 

C/(r) = Uo{r) + AC/,ct.(r) + AC/,efl.(r) + 0{T-^) (2.16) 

with 



AC/,(r) 
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ReAPl,(r), (2.17) 



i — ret.,refl. The relative corrections due to retardation 
and reflection read 
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Note that ReAPtj(r) is an even function of w, as fol- 
lows directly from definition (j2.5p together with the 
Schwarz reflection principle G(r,r';— w) = G*(r,r';a;). 
As a consequence, the leading temperature corrections 
in the high-temperature limit, being proportional to 
ReAP(^(r)/a;, change sign when comparing downward 
and upward transitions. 



III. CASIMIR-POLDER POTENTIAL NEAR A 
SPHERE 

As depicted in Fig-IH we consider a particle at distance 
r from the center of a sphere of radius R and permittivity 
£ = e{(jj). The dyadic Green's function leads to [37| 

00 

1=1 
+ r™{^x) [lil + l)h\'\xf + h'i''^\xf] } (3.1) 

where h\ ' (x) is the spherical Hankel function of the first 

kind, hi (x) [and for future reference, J;(2;)] is shorthand 
for 

h<i'^'ix)^[xh'i'\x)y; j\ix)^[xMx)]' (3.2) 

[ji{x): spherical Bessel function of the first kind]. For 
convenience we are using the dimensionless distance and 



size parameters x = ruj/c and Q < (j) — R/r < 1, recall 
Eqs. pn|) and ([TTB)) . The reflection coefficients for TE 
and TM-polarized waves read 
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(3.3a) 



(3.3b) 



In the following we assume both the particle-center 
separation and the sphere size to be non-retarded, (fix < 
a; ^ 1. In addition, we perform the perfect-conductor 
limit |e| :^ 1. To leading order, the two limits commute: 
Taking the perfect-conductor limit first, the reflection co- 
efficients reduce to 



„TE 



(0x) — > r; ((px) 



,TM(^^)e^^TM.PC(^^) 



'h['\.pxy 

hl'^'ict^x) 



(3.4) 



(3.5) 



see the asymptotes (jAll) and (jA2p in App. [X] Using the 
expansions (IA3[) and (|A4[) . they further simplify to 
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/ (2/ + l)!!(2;-l)!! 
(3.7) 



in the non-retarded limit. Here, (2/-fl)!! = 1-3 ••• (2^+1). 
In contrast, when taking the non-retarded limit first, 
one finds 
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and subsequently 
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(3.11) 



While the TM-coefficicnt takes the same form regardless 
of the order of the limits, the TE-coefficient yields dif- 
ferent results, depending on which of the two limits is 
performed first. 

However, the Green tensor in the non-retarded limit is 
dominated by r™ . Substituting the results for the reflec- 
tion coefficients into p.ip , using the approximation (|A4[) 



from App. |X] and retaining only the leading order in x, 
one finds 

-. oo 



1 (/)^(6-30^-H(/)^) 
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With this result, the temperature-invariant CP poten- 
tial (pn^ reads 
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(3.13) 
-^ 1. 



(3.14) 

This is in agreement with the zero-temperature potential 
as found in Ref. [38| for a perfectly conducting sphere in 
the non-retarded limit on the basis of image-charge tech- 
niques. As discussed in Ref. [43|, the atom-sphere geome- 
try is a particular example of a two-parameter geometry, 
conveniently described by a scaling function /(0). 

The accuracy of the simplest, temperature- 
independent approximation, Eq. p.l3p . is demonstrated 
numerically in Fig. [2] (dashed lines) for a ground-state 
two-level particle outside a gold sphere. The permittivity 
of the sphere has been described by a Drude model 
e{uj) = 1 — ujp/[uj{uj + ij)] with parameters ujp = 9eV 
and 7 — 35meV. For comparison, the same situation 
but with a smaller sphere is shown in Fig. |3l We 
see that Eq. p.l3p yields a very good approximation. 
The exact potential is slightly smaller by at most 5% 
for rcj/c = 0.1. Recall that the leading correction 
in the high-temperature limit has opposite signs for 
ground-state and excited atoms. For an excited two-level 
atom, we would hence find that the exact potential 
is slightly larger than its approximation p.l3p . In 
the following, we derive analytical expressions for the 
leading temperature-dependent corrections to Eq. p.l3p 
providing an even much improved approximation at 
higher temperatures. 



Correction from retardation 



When considering the full thermal CP potential ()2.2ap 
and (j2.2bp using the Green's function of Eq. (13.11) . the 
non-retarded and perfect-conductor limits do not com- 
mute. We have 0a; < x ^ 1 and |e| 3> 1, which is com- 
patible with both large and small values of |-\/e|0a;. For 
a metal sphere at typical experimental distances of order 
micrometers and x ~ 0.01 — 0.001, we have Im-ye0a; >> 1, 
meaning that the perfect-conductor limit has to be ap- 
plied first. The opposite limit |v^|0a: <C 1 may be re- 
alised for dielectrics whose permittivity tends to some 
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FIG. 2: Numerical comparison of the exact CP potential U{T) 
of the contribution from a transition of energy hck for a par- 
ticle outside a gold sphere with f/approx being respectively the 
T-independent result Uo for a perfect conductor in the non- 
retarded limit (dashed line) and the approximate expression 
Eq. (|1.5p including the linear temperature correction (solid 
line). Parameter values; R = ir = lO/xm. 



moderate electrostatic value, in which case the nonre- 
tarded limit would have to be performed first. We briefly 
consider this case in Appendix [BJ 

The correction from retardation effects is found by us- 
ing the perfect conductor values p.4p and (I3.5P of the 
reflection coefficients and expanding rjj*-'(r) as given by 
Eq. p.l|) in powers of a;, the leading correction term being 
of order x^ . We obtain such quadratic corrections from 
three sources: (A) from the TM-mode reflection coeffi- 
cient Tj ' ; (B) from TM-mode propagators /ij (x)^ 

and h\ (x)"^; and (C) from the leading-order TE-mode 
contribution. The technical details of the small x expan- 
sions of the different cases are found in App lA ll 
As shown therein, correction (A) takes the form 
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with r, 
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being given by Eq. (j3.7p . The correction 



(B) is found to be 
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FIG. 3; Same as figure (2] but with a smaller gold sphere: 
R = ^r — Ifim. 



Finally, the correction (C) from the TE-mode takes the 
simple form 



TE,PC/, x,(l). ^2 <^^'"^^ 


+ ... 


(3.17) 
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Substituting these corrections into Eq. (]3.il). 


one finds 
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The sum may be carried out in closed form by splitting 
the expressions into partial fractions. One finds 
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Substituting this result into Eq. (j2.17p . we obtain the 
retardation correction 

with scaling function 

5rct.((/') 

= 3(1 + 3(/)^)artanh0 - 0(3 - 0^) + 20^ log(l - (j)'^) 

-.1'^' , ^°^^-0' (3.21) 

l-41og(l-0) for 0^1. ^ ^ 



Its relative contribution (|2.18l) is given by 



Uoir) 



tvjj 



gret.(^) 
/(0) ■ 



(3.22) 



B. Correction from imperfect reflection 

The leading order corrections to the ideal reflection 
coefficients are calculated in App. lA2l and have the form 



TE,PC/ 



rl^c^x) =r,\^''-^i^x) 



TM,PC/ 



r/-(<^a;)=r,7'--(</.x) 



i{2l + l) 
y/e4>x 
i(hx 21 + 1 



Fe l{l + l) 



(3.23a) 
(3.23b) 



for X <^ 1, with Tj Q ' and r^ g ' given by Eqs. (I3.6p 
and p.7p . Substituting these results into Eq. p.ip and 
using the leading-order small argument expansions in 
Eqs. ([XTC]) and (jXTTI) . we find 



V 1 — 1 



ix<p 



Airr'^y/e 



1=1 

2 aa4: 



21 + 1 



,21 



3 + 7r - A(j) 

(1-02)2 



- log(l - cf") 



(3.24) 



Inserting this into Eq. ()2.17p . the refiectivity correction 
is found to be 



A[/rcfl.(r) 



\dVx 



.Rc(^ 



kv.T 1 
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with scaling function 

'3 + 702 _ 4^4 



hw 



■^j 5rcfl.(0) 

(3.25) 



'602 



(1 



(1-02)2 

for -> 0, 
for — > 1. 



l0g(l - 02) 



(3.26) 



Its relative contribution p.lSp reads 
AC/refl.(r) _ fkBT 1 



Uo{r) 



^iw 2J'''^°U/^j /(0) 



(3.27) 



C. Discussion and comparison 



Combining the T-invariant result (j2.12p with the re- 
tardation and reflectivity corrections (I3.20p and p.25p . 
we obtain the weakly temperature-dependent CP poten- 
tial (|1.5p . as stated in the introduction. The quality of 
this analytic high-temperature result is demonstrated in 
Figs. [5] and [3] (solid lines), where we compare it with the 
result of an exact numerical calculation for the contribu- 
tion Unk [see Eq. (II. 2p ] from an upward internal energy 



100 



gW/m) 



0.01 



0.001 



Small 
curvature 



Large 
curvature 




0.0 



FIG. 4: Ratios of the scaling functions (?rct.(0)//(0) (solid 
line) and grefl. (0)//(0) (dashed line) that govern the impact of 
retardation and finite refiectivity on the thermal CP potential. 



transition of a particle outside a gold sphere. One sees 
that the analytic result is an excellent approximation for 
temperatures T > 200 K. Notice in particular that for x 
of the order 0.1 and lower, Eq. (|1.5p is an excellent ap- 
proximation (better than 1%) at T « 300K, where most 
experiments are performed. 

While the T-invariant flrst term in Eq. (jl.5|) is quite 
adequate for very small retardation values (such as a; = 
0.001, typical of Rydberg atoms), the full expression is 
much better as x increases to about 0.1. This could be the 
case for certain cold polar molecules. For example, with 
LiH molecules [20], x = 0.1 corresponds to r = 15/im 
for the dominant, rotational transition which is not an 
atypical situation. In this case the T-independent term 
is off by almost 10% at 300K with the 10 micron sphere, 
but is still better than 1% when the correction is in- 
cluded. In general, the T-invariant first term in Eq. (jl.Sp 
becomes a worse approximation at higher temperatures, 
with the error increasing without bounds as the envi- 
ronment temperature rises. In contradistinction, our ap- 
proximation (|1.5p becomes better and better at high tem- 
peratures and its error remains bounded throughout. 



According to Eqs. p.22p and (|3.27p . the relative con- 
tributions from retardation and finite refiectivity are gov- 
erned by the ratios gi{4>)/ f{<p) of the scaling functions as 
given by Eqs. ([XTi)) . (P%t and (15:^ . These ratios are 
depicted in Fig. 01 The figure shows that both contribu- 
tions strongly depend on the curvature of the sphere as 
parametrised by 0. The retardation contribution takes a 
value 1/3 for a strongly curved sphere and stays approxi- 
mately constant for < 0.5. In the limit of a fiat surface, 
it rapidly falls off as — 4(1 — 0)^ log(l — 0). The refiec- 
tivity correction grows as 1/0 in the limit of a strongly 
curved sphere and falls off gently as 6(1 — 0) in the limit 
of a flat surface. 

The ratio grofl.(0)//(0) is greater than grot. (0)//(0) by 
at least an order of magnitude for all curvatures. One has 
to bear in mind, however, that the reflectivity correction 
carries an additional factor Ke{i/y/e)/x ^ 1. For a metal 
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described by a Drude model with w, 7 ^ top, one finds 



J y/u}^/2/ujp for oj ^ 7, 



\^/ 



wp 



for oj 3> 7. 



(3.28) 



Depending on the actual value of Re{i/y/e)/x ^ 1 for a 
given molecule and material, either one of the reflectivity 
or retardation may dominate for given curvatures. How- 
ever, the asymptotic behaviour observed in Fig. |4] shows 
that the reflectivity correction will always dominate in 
the limits of small or large curvature. 



IV. SUMMARY 

We have studied the temperature-dependent CP po- 
tential of a particle strongly out of thermal equilibrium 
near a metal sphere, and have derived a simple approx- 
imate expression in closed form for the interaction po- 
tential. The approximation is valid at the 1% or bet- 
ter for Rydberg atoms and cold polar molecules at all 
temperatures for typical experimental length scales. We 
have assumed both the particle-sphere distance and the 
particle's transition frequency to be small enough so 
that retardation and imperfect reflectivity present only 
small perturbations to the temperature-independent re- 
sult. Such is typically the case in experimental set-ups in 
which cold polar molecules or Rydberg atoms are used, 
whose interaction potential is dominated by long wave- 
length transitions for which the non-retarded regime ex- 
tends to tens and hundreds of micrometers, respectively. 

In recent publications it has been shown that the 
Casimir-Polder potential acting on a particle prepared in 
an eigenstate at a non-retarded distance from a macro- 
scopic body can be virtually independent of the sur- 
rounding temperature. This is the case for the geom- 
etry considered, and the error made in approximating 
the interaction as independent of temperature from zero 
to room temperature is only a few percent for suffi- 
ciently non-retarded interaction. The small temperature- 
dependent corrections to the potential of a metal body 
have been identified to stem from retardation and imper- 
fect reflectivity and we have analysed these separately 
and discussed the relative importance of each with re- 
spect to the other. Our results show that reflectivity is 
the dominant correction for very large or small curva- 
tures, while intermediate curvatures may be governed by 
either retardation and reflectivity corrections, depending 
on particle and material. 

The perturbative method employed in this investiga- 
tion is equally well suited for the study of more com- 
plicated geometries. Again, it promises physical insights 
that are hard or even impossible to gain by numerical 
means. 

We thank Ho Trung Dung for discussions. This work 
was supported by the UK Engineering and Physical Sci- 



ences Research Council. Support from the European Sci- 
ence Foundation (ESF) within the activity 'New Trends 
and Applications of the Casimir Effect' is gratefully ac- 
knowledged. 



Appendix A: Limits and their leading corrections 

To calculate the CP potential in the perfect-conductor 
and non-retarded limits, we need to approximate the 
spherical Bessel functions for small and large arguments. 
Using the asymptotes from §10 of Ref. [43] 

. , . sin(a; - ln/2) 
ji{x)k. — i '—L forx>l; (Al) 



hf\x) 



(_i)'+le'(2:-«'r/2) 



for a; ^ 1, 



(A2) 



we easily find the perfect-conductor limits p.4[) and (13.51) . 
The non-retarded limits dSH), dSj]), dSH) and ((3:9)) 
can be found by using the expansions [43i] 



ji{x) 



ftW(x) 



(2; + !)!! 
1(2/ -1)!! 



1 



2(2/ -F 3) 
^ 2(2/ -1) 



+ 



(A3) 
(A4) 
Note 



for X < 1 with (2/ + 1)!! = 1 • 3 • 5 • • • (2/ -I- : 

that the next-to-leading order expansion is only needed 

for (|3.8p . where the leading-order term vanishes. 



1. Retardation corrections 

To calculate the retardation correction to the perfect- 
conductor TM- mode reflection coefficients p.5|) . we use 
the expansions (IA3p and (jA4p . which together with the 
definitions (13.21) lead to 



~3'i{x) 



;(!)' 



{x) 



(/ + l)x' 
'(2/ + 1)!! 
i/(2/-l)!! 



1- 

H 



l + i 



2 (2/ + 3)(/-Hl) 
a;2 1-2 



2 l{2l~ 1) 



(A5a) 
(A5b) 



for a; ^ 1. This immediately yields Eq. (I3.15p . 

The retardation correction from the propagator factors 
are found from expansions 



h^^\xf 



[(2/-l)!!]2 



„2/+2 



1 



;(i)' 



[xf 



/^[(2/-l)!!]^ 
^21+2 



2/- 1 

:c'(/-2) 
^ /(2/-1) 



(A6a) 
(A6b) 



which follow from Eqs. (IA4[) and (jA5bp . Combining these 
results with r, ^ ' as given by Eq. p.7p . we arrive at 



';,o 



Eq. (IXTCl) . 

Expansion (|A6aP moreover combines with r^^' as 
given by Eq. (P?B| to resuh in Eq. (|XT7)) . 



2. Finite reflectivity correction 

In order to expand the reflection coefiicients p.3p in 
powers of e~^, we rewrite them as 



„TM/ 



hx) =r 



TM.PC/ , \ 1 ^ AJ/e 






where 



^^jjWf^. j^ji{4>x) 



ji{^/e(j>x) 



fiifl^x) 



H ^ 






With the assumption Iml^/ej^x ^ 1, the asymp- 
tote (PTj) leads to 



j;(Ve0a;) 



:^\/e0a;cot(\/£0a:; — ^) « ~iy/e(j)x. 



jiiVs(t>x) 
Using the asymptotes (jASp and (jA4l) . we further have 

1 



/.P)(0x) 



for a; <Si 1. Combining these results and retaining 
only the next-to-leading order in x, one easily obtains 
Eqs. dSSl). 



Eq (jBip is the counterpart of (|1.5p for the case of a di- 
electric sphere. 



As for the case of a metal sphere, the approxima- 
tion is good to at least the 1% level at all tempera- 
tures when X < 0.01, and becomes significantly better 
ikr=0.001 




50 100 
Temperature (K) 



500 1 000 



FIG. 5: Numerical comparison of the exact CP potential U{T) 
of a two-level model particle (transition energy hck) outside a 
dielectric sphere with the approximation Eq. (IBlf) with (solid 
line) and without (dashed lines) the correction term in the 
curly brackets. Parameter values: R= ^r = lO^m and e = 6. 



Appendix B: Casimir-Polder expression for a 
dielectric sphere 



An expression for the CP potential on a particle near 
a dielectric sphere can be readily derived using the same 
methods as elsewhere in this article. Although it may be 
written in closed form in terms of hypergeometric func- 
tions, this hardly constitutes a simplification, and we will 
only give the expression as an infinite sum. The sphere's 
dielectric constant e is now no longer assumed large, so 
that Iv^la; ^ 1 is assumed. The TE contribution is now 
of order x"* and can be neglected. Since we are con- 
sidering frequencies Wkn which are typically small on an 
optical scale, we assume e(wfc„) = £(0) — e. The result is 



rrdicl. 



|dfe„|2(e-l)^?(/ + l)(2l + l)</>2'+i 



247reor'' 
{l + 2a;2( 



1=1 
k^T 1 



el + 1 + 1 



hwk 



21 + 1 



\2l + l) e{l-2) + l + l 



{2l-l){2l + Z) el + 1 + 1 



}• (Bl) 



for T > hhj/kB (about UK for x 
bers in figure [S|) . 



0.01 with the num- 



When assuming as we have that e does not vary appre- 
ciably between frequencies and ojfcn, the leading order 
correction is of order x^. It is noteworthy that the T- 
independent term alone (without the correction in the 
curly braces of Eq. (|B1I) ) is about as good an approxima- 
tion as the corresponding expression is for a metal sphere. 
This conclusion would not hold, however, if the dielectric 
material has resonances at frequencies in the order of w, 
which might lie in the microwave or far infrared regime, 
in which case AF^j may no longer be small compared to 
Fq (see Section HI]) . 



[1] H. B. G. Casimir and D. Polder, Phys. Rev. 73, 360 
(1948). 



[2] T. F. Gallagher, Rydberg Atoms (Cambridge University 



10 



Press, 1994). 
[3] H. Kiibler, J. P. Shaffer, T. Baluktsian, R. Low, and 

T. Pfau, Nature Photonics 4, 112 (2010). 
[4] A. Tauschinsky, R. M. T. Thijssen, S. Whitlock, H. B. 

van Linden van den Heuvell, and R. J. C. Spreeuw, Phys. 

Rev. A 81, 063411 (2010). 
[5] A. S. S0rensen, C. H. van der Wal, L. L Childress, and 

M. D. Lukin, Phys. Rev. Lett 92, 063601 (2004). 
[6] P. Hyafil, J. Mozley, A. Perrin, J. Tailleur, G. Nogues, 

M. Brune, J. M. Raimond, and S. Haroche, Phys. Rev. 

Lett 93, 103001 (2004). 
[7] D. Petrosyan, G. Bensky, G. Kurizki, I. Mazets, J. Ma- 

jer, and J. Schmiedmayer, Phys. Rev. A 79, 040304(R) 

(2009). 
[8] S. Y. T. van de Meerakker, H. L. Bethlem, and G. Meijer, 

Nature Physics 4, 595 (2008). 
[9] L. D. Carr, D. DeMille, R. V. Krems, and J. Ye, New J. 

Phys. 11, 055049 (2009). 
[10] M. T. Bell and T. P. Softley, Mol. Physics 107, 99 (2009). 
[11] S. A. Meek, H. Conrad, and G. Meijer, Science 324, 1699 

(2009). 
[12] A. Ashkin, Phys. Rev. Lett. 24, 156 (1970). 
[13] J. Chen, J. Ng, Z. Lin, and C. T. Chan, Nature Photonics 

5, 531 (2011). 
[14] S. Kawata and T. Sugiura, Opt. Lett 17, 772 (1992). 
[15] S. M. Spillane, T. J. Kippenberg, and K. J. Vahala, Na- 
ture 415, 621 (2002). 
[16] H. Haakh, F. Intravaia, C. Henkel, S. Spagnolo, R. Pas- 

sante, B. Power, and F. Sols, Phys. Rev. A 80, 062905 

(2009). 
[17] S.Y. Buhmann, M.R. Tarbutt, S. Scheel, and E.A. Hinds, 

Phys. Rev. A 78, 052901 (2008). 
[18] J. A. Crosse, S. A. EUingscn, K. Clements, S. Y. Buh- 
mann, and S. Scheel, Phys. Rev. A 82, 010901(R) (2010); 

Erratum ibid. 82, 029902(E) (2010). 
[19] S. Y. Buhmann and S. Scheel, Phys. Rev. Lett. 100, 

253201 (2008). 
[20] S. A. EUingsen, S. Y. Buhmarm, and S. Scheel, Phys. 

Rev. A 79, 052903 (2009). 
[21] S. A. EUingsen, S. Y. Buhmann, and S. Scheel, Phys. 

Rev. Lett. 104, 223003 (2010). 
[22] S. A. EUingsen, S.Y. Buhmann, and S. Scheel, Phys. Rev. 

A 84, 060501(R) (2011). 
[23] T. Nakajima, P. Lambropoulos and H. Walther, Phys. 

Rev. A 56, 5100 (1997). 



[24] S.-T. Wu and C. Eberlein, Proc. R. Soc. Lond. Ser. A 
456, 1931 (2000). 

[25] M.-P. Gorza and M. Ducloy, Eur. Phys. J. D 40, 343 
(2006). 

[26] M. Antezza, L. P. Pitaevskii, S. Stringari, and V. B. Sve- 
tovoy, Phys. Rev. A 77, 022901 (2008). 

[27] Y. Sherkunov, Phys. Rev. A 79, 032101 (2009). 

[28] S. Y. Buhmann and D.-G. Welsch, Prog. Quantum Elec- 
tron. 31, 51 (2007). 

[29] S. Scheel and S. Y. Buhmarm, Acta Phys. Slov. 58, 675 
(2008). 

[30] S. A. EUingsen, Y. Sherkunov, S.Y. Buhmann, and 
S. Scheel, in Proceedings of the Ninth Conference on 
Quantum Field Theory Under the Influence of Exter- 
nal Conditions (QFEXT09) edited by M. Bordag and 
K. A. MiUon (World Scientific, 2010), p. 168, Preprint: 
quant -ph/0910 . 5608. 

[31] A. Weber and H. Gies, Phys. Rev. Lett. 105, 040403 
(2010). 

[32] T. L. FerreU and R. H. Ritchie, Phys. Rev. A 21, 1305 
(1980). 

[33] A. M. Marvin and F. Toigo, Phys. Rev. A 25, 782 (1982); 
25, 803 (1982). 

[34] C. Girard, S. Maghezzi, and F. Hache, J. Chem. Phys. 
91, 5509 (1989). 

[35] W. Jhe and J. W. Kim, Phys. Rev. A 51, 1150 (1995). 

[36] V. V. Klimov, M. Ducloy, and V.S. Letokhov, J. Mod. 
Opt. 43, 2251 (1996). 

[37] S. Y. Buhmann, H. T. Dung, and D.-G. Welsch, J. Opt. 
B: Quantum Semiclass. Opt. 6, S127 (2004); Erratum: 
J. Phys. B: At. Mol. Opt. Phys. 39, 3145 (2006). 

[38] M. M. Taddei, T. N. C. Mendes and C. Farina, Eur. J. 
Phys. 31, 89 (2010). 

[39] A. Sambale, S. Y. Buhmann, and S. Scheel, Phys. Rev. 
A 81, 012509 (2010). 

[40] M. Babiker and G. Barton, J. Phys. A: Math. Gen. 9, 
129 (1976). 

[41] L Brevik, S. A. EUingsen and K. A. Milton, New J. Phys. 
8, 236 (2006). 

[42] S. Y. Buhmann, S. Scheel and J. Babington, Phys. Rev. 
Lett. 104, 070404 (2010). 

[43] M. Abramowitz and L A. Stegun, Handbook of Mathe- 
matical Functions (Dover, New York, 1964). 



